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Abstract 

A new solution to the star-triangle relation is given, for an Ising type model of interacting 
spins containing integer and real valued components. Boltzmann weights of the model are given 
in terms of the lens elliptic gamma function, and are based on Yamazaki’s recently obtained 
solution of the star-star relation. The star-triangle relation given here, implies Seiberg duality 
for the 4 —d Af = 1 Si x S^jZr index of the SU{2) quiver gauge theory, and the corresponding 
two component spin case of the star-star relation of Yamazaki. A proof of the star-triangle 
relation is given, resulting in a new elliptic hypergeometric summation/integration identity. The 
star-triangle relation in this paper contains the master solution of Bazhanov and Sergeev as a 
special case. Two other limiting cases are considered one of which gives a new star-triangle 
relation in terms of ratios of infinite g-products, while the other case gives a new way of deriving 
a star-triangle relation that was previously obtained by the author. 
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1 Introduction 


The star-triangle relation is a distinguished form of the Yang-Baxter equation for Ising-type models 
on two-dimensional lattices. In these models the fluctuating variables, or “spins”, are assigned to 
lattice sites, and two spins interact only if they are connected by an edge of the lattice. Remarkably, 
many physically interesting models in this class can be solved exactly, for instance, the 2—d Ising [T], 
and chiral Potts mm models, and some others UHlo] (see also mm for a review of other known 
cases). The star-triangle relation plays the role of the integrability condition for these models. 

Recently Bazhanov and Sergeev (BS) obtained an important “master” solution [lOj of the star- 
triangle relation, which contained all previously known solutions of this relation as particular cases, 
and provides interesting new examples. The above master solution is expressed in terms of the 
elliptic gamma function, which contains two arbitrary free parameters p and q, that play the role of 
elliptic nomes. The spin variables for the corresponding statistical mechanical model take continuous 
real values on the circle. 

Considered as a mathematical identity the BS master solution is identical to the elliptic beta 
integral of Spiridonov m- The latter discovery was central to the modern development of the theory 
of elliptic hyper geometric functions m, and some recent works further highlight that some of these 
identities are connected to the integrability of lattice models of statistical mechanics. Some examples 
include an extension of the BS master solution to the case of multi-component spins [15(116]. and 
remarkable correspondences to Seiberg duality in supersymmetric gauge theories UZHlQ]. Recently 
Yamazaki introduced a new integrable model [20], with Boltzmann weights satisfying the star-star 
relation, by using the property that the latter relation is equivalent to a particular Seiberg duality for 
the 4—d J\f = 1 lens index for a class of SU(N) quiver gauge theories. This star-star relation is rather 
general and contains its variant for the master solution m and its multi-spin generalisation [151116] 
as particular cases. 

In Section [2] it is shown that the Boltzmann weights for the model with two-component spins 
introduced by Yamazaki, also satisfy a star-triangle relation. A proof for the star-triangle relation 
is given in AppendixjA] which also verifies the corresponding two component spin star-star relation, 
since the former relation implies the latter (but the reverse is not true). The actual proof given in 
AppendixjAjis for an identity more general than the star-triangle relation, resulting in a new elliptic 
hypergeometric summation/integration identity for the lens elliptic gamma function, that contains 
six complex and six integer variables. This identity contains Spiridonov’s celebrated elliptic beta 
integral as a particular case. Two limiting cases of the star-triangle relation are considered, resulting 
in one new solution to the star-triangle relation with Boltzmann weights given in terms of infinite 
q-products, and another new solution with Boltzmann weights given in terms of the Euler gamma 
function, that has recently been obtained by the author [2T]. Possible relations to existing integrals 
in the literature are discussed. 

1.1 Solvable square lattice model 

All models may be considered here on the square lattice made up of of N sites. Spin variables 

aj = {xj,mj ), Xj € M., ruj £ Z , j = 1,2,..., N, (1) 

are assigned to each site of the lattice, where Xj takes real values, and rrij takes integer values. Two 
spins interact only if they are connected by an edge of the lattice. The interactions are represented 
by the Boltzmann weights Waio'i,aj), and Waio'i,crj), associated to horizontal and vertical edges 
respectively, where crj and aj are the spins located at the end of the edge, as shown in Figure [H 
Here two Boltzmann weights are distinguished by crossing of dashed rapidity lines, a property which 
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Figure 1: Horizontal (left) and vertical edges (right), and their Boltzmann weights. Here “rapidity lines” 
(dashed arrows) are shown to distinguish the two types of weights, W and W. 


allows one to also consider the model on more general “Z-invariant” lattices [22]. The two edge 
Boltzmann weights, depend on the additive spectral variable a, and are related to each other by 
the crossing symmetry property Wa(<7i,(Tj) = Wrj-aio'i,crj). The “crossing parameter”, rj, is model 
dependent, and the regime 0 < a < ?? is a physical regime of the model, where the Boltzmann weights 
are positive and real-valued. For all models considered in this paper, the Boltzmann weights are 
also spin reflection symmetric, such that Wa(cTi,fTj) = Wa(fTj, Uj). 

To each spin aj in the lattice one also associates the single-spin weights S{aj), which are inde¬ 
pendent of the spectral variable a. The partition function of the model is then defined as a product 
of all Boltzmann weights, with a integral (sum) over all internal continuous (discrete) spins, while 
boundary spins are kept fixed, 

^=E/n FVq,( fTj, fJjf) Wrj— Q, (fTfc , IT;) ‘S(o'n) dXn ■ (2) 

(b) m 

The hrst product is taken over all horizontal edges (ij), the second over all vertical edges (kl) and the 
third product over all internal sites of the lattice. The goal of statistical mechanics is to evaluate ([2|) 
in the thermodynamic limit, when the number of sites in the lattice goes to inhnity. This evaluation 
is possible if the Boltzmann weights satisfy the star-triangle equation. For the models given here 
this relation reads 

^ / dxo5(cro)Wr,-ai(fTi,0-o)W^_„.((Tj,Cro)Wr,-a;,(fTo,(Tfc) 

mo (^) 

= 7Z{ai,aj,ak) Wq,. ((Tj,ak)yVaj {(Ti,ak)Wai, {crj,(Ti) , 

where the three spectral parameters a;, aj, ak satisfy the constraint ai + Oj + ak = rj and the factor 
'lZ{ai, aj,ak) is independent of the spins Uj, crj,ak- The integral and sum are evaluated over a given 
set of continuous and discrete values respectively. There exists also a second star-triangle relation 
obtained by exchanging the order of spins appearing in the Boltzmann weights in ([3|), however 
for models with symmetric Boltzmann weights, Waio'i,o'j) = Wo((Tj,Uj), the two expressions are 
equivalent. The star-triangle relation ([3|) implies that the row-to-row transfer matrices of the lattice 
model commute [23]. 

For all models considered here, the normalisation of the Boltzmann weights is chosen such 
that the spin independent factor Tl{ai,aj,ak) in ([3]) is equal to one. This result is based on a 
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factorisation for TZ{ai,aj,ak) |12j . which holds for all of the above mentioned solutions of the 
star-triangle relation, particularly for the new solutions appearing in the following section. 

For this special normalisation, the Boltzmann weights of the model also satisfy the following 
boundary conditions 

VV’a(c’'i, crj)|^_g = 1, + —Xj) } (4) 

where d{x), and 6m,n, are respectively Dirac and Kronecker delta functions. The exact form of the 
second boundary condition differs slightly depending on the symmetries of the Boltzmann weights, 
and explicit expressions for the boundary conditions in each case will be given for the three different 
models obtained in the next section. 

From the boundary conditions Q, and star-triangle relation ([3]), one obtains the following 
inversion relations 


mo 


yVa{(7i,Crj)W-a{(Ti,(Tj) = 1 , 

dxoS{ao)WrJ-a{<7i,(Xo)y^1^+a{cro, CTj) = 


((f(Xj-|- Xj') 5mi,—mj T 6(yXi Xj) 6mi,mj') ■ 


25 (a, 


The above relations, ([3]) and ([5]), allow one to show that in the thermodynamic limit as the number 
of lattice sites goes to infinity N ^ oo, the bulk free energy of the model vanishes 

lim Ai“^log.E = 0. (6) 

N—^OO 

A derivation of this result requires some extensions of the standard inversion relation method 
|24H26j . Here the boundary spins are assumed to be kept finite in the limit N ^ oo, and there is 
an analyticity assumption for the free energy of the model in the physical regime. The result ([6]) is 
purely a consequence of the special choice of normalisation for the Boltzmann weights [HHlo]. 


2 New discrete and continuous spin solutions to the star-triangle 
relation 

In this section the Boltzmann weights are defined that give a new solution of the star-triangle 
relation (l3])0 This star-triangle relation and the corresponding proof given in Appendix are the 
main result of the paper. 

Recall the definition of the spin 

aj = {xj,mj), Xj G M , rrij G h . (7) 

Restrict the continuous real valued component Xj, and the discrete integer valued component mj, 
to take values 

0 < Xj < vr, ruj = 0,1,..., [r/2j , (8) 

for some positive integer parameter r = 1,2,..., where [ J is the floor function. Dehne also the 

elliptic nomes p, q, and crossing parameter r] as 

p _ gTTio-, q _ gTTir, ^ _ _^j(^ _ 1 _ ^)/2 , Imcr, lmr>0. (9) 

^ These Boltzmann weights correspond to the 2-component spin version of Yamazaki’s star-star relation, which 

was obtained by identifying the correspondence of this relation with Seiberg duality of the lens index for a, Af — 1 
supersymmetric quiver gauge theory m- 
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Note that a physical regime where Boltzmann weights are real and positive valued can be found for 
p = q*. Define the elliptic gamma function as 


( 10 ) 


I _ q2\z p2jr+l q2fc+l 


n 1 _ e-2i2 p2i+i q2fc+i • 


j,k=0 

In terms of the elliptic gamma function, the so-called lens elliptic gamma function is defined as 
^r,miz) = (r /2 - lmjr)'^o-;pq, p^) $(2 - (r /2 - Iml,.) vrr; pq, q'’) 


n^T 


1 _ p-2Mr ('pq^ 2 j-i-i |'pr^2fc-i-2 2 — q2Mr ('pq) 2 j-i-i ('qr'^2fc 


( 11 ) 


j,k=0 


— p2Mr (pq)2j-i-i (p»')2fc I _ e- 2 i 2 : q-2[m]r (pq)2i+i (q»-)2fc-i-2 ’ 


where |m]]r € {0,1,... , r — 1} denotes m modulus r. From ([ 8 ]), note that when r = 1, then rrij = 0, 
and the lens elliptic gamma function reduces to the usual elliptic gamma-function (| 10 p 

^>i,o(z) = ^>(z;p,q). ( 12 ) 

The lens elliptic gamma function (|lip . satisfies the following periodicity and inversion relations 

1 


^r,m{,z) — ^r,m(^ T 

Now define the edge Boltzmann weight as 




— m( z') . 


(13) 




K{a) 


*^r,mi—mj (,Xi Xj icxj ^r,mi+mjixi Xj \u) 


where |m]± := |m]]r|—mjr- The spectral parameter a is taken to lie in the domain 0 < a < 77 , 
where r] is real. For p = q*, this is a physical regime of the model, where the Boltzmann weights 
(fill) take real, positive values. 

The normalisation factor K{a) is given by 


K{a) = exp 


:.4Q:n 


((pq)™ - (pq)-™) 


^ n((pq)2^ - (pq)-2^)(p™ - p-™)(q™ - q-^*^) 


and satishes the pair of functional equations 

K{r] — a) 


K,{a) 


= <kj,,o(i(^ “ 2 a)), K{a)K{—a) = 1 . 


(15) 


(16) 


For r = 1 this reduces to the normalisation of the Boltzmann weights for the BS master solution HUj. 

Note that the functional equations (1161) arise when solving for the free energy of the model, 
hm 7 v_>.oo log 2, using the inversion relation method |241426] . The solution of these functional 
equations with appropriate analyticity properties, (IISl), is included in the normalisation of the 
Boltzmann weights (jl4p . and the result ([6]) follows (SHUj. 

Next define the single-spin Boltzmann weight aqj 

5(a,) = ^ (p2’'; p 2 ’')oo(q'";q'’')oce 2 '^Pm,l±A - \r,) $,, 2 „.,( 2 x, - \r ,), 

= — 7?4(2xj -h (r/2 - |2mi]r)7ro-1 p'’) 'd 4 ( 2 xi - (r/2 - pmj]]r)7rr | q'"), 


^This differs from Yamazaki’s single-spin weight (S^fs) in his notation [^) by dropping a constant singular factor 
that appears to be incorrect, at least from the point of view of the statistical mechanical model (in his notation this 


factor is (Fr.ofl;p,g)) 


-(JV-l)x 
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where 



is a Jacobi theta function 


mi = 0 or |r - niijr , 
otherwise, 


Mz I P) = (P^; p')oo n (1 - (1 - , 

n=l 


and (x;q)oo = nr=o (1 — xq^) is the q-Pochhammer symbol. 
The Boltzmann weights (jl4ll are reflection symmetric 


(18) 


(19) 


Wa((Ti, fJj) = Wa{(rj,(Ti) . (20) 

The Boltzmann weights have an obvious vr-periodic symmetry in the continuous spin variable, and 
they also are invariant under the spin transformation Xi —>■ —Xi, nii —>■ r — rrii, Accordingly the 
discrete spins are restricted to values 0,1,... , [r/2j, and the e* factor was introduced in (fTSl) to 
account for this. 

The Boltzmann weights (jl4p satisfy the following boundary conditions analogous to ([1|) 


y^a{^h^j)\n,=n = 


bV»7—a(^i) ^j)|q_^q S[(t ) ^[rni+mjlrfl (5(sin(xj Xj)) mj]r,o) ) 


( 21 ) 


where S{ai) ^ 0. The r = 1 case of these relations were previously obtained for the BS master 
solution m, and in connection with biorthogonality of elliptic hyper geometric functions |28j . 

The Boltzmann weights (jl4p . and (I17p . satisfy the star-triangle relatiorU 


/ dxoS{ao)W^-aii(ri,ao)Wr,-aj{(^j,(^o)yVv-aki^k,(^o) = Vya.(<Tj,CTfc)Wa,(<Ti,afc)>Vaj^(fTj,CJi) , 

mo=0 

( 22 ) 

with the spectral parameters satisfying rj = ai + aj + ak- For r = 1 this reduces to the master 
solution of the star-triangle relation [10] . 

The star-triangle relation (I22p is a particular case of a new elliptic hypergeometric summa¬ 
tion/integration identity given in Appendix lAl 


2.1 Limit: r —)■ cxd 

The r —>■ oo limit of (I22p is formally fairly straightforward due to the simple asymptotics of the 
lens elliptic gamma function. Consider the same elliptic nomes p, q from the previous section in Q. 
Dehne the function Q as the r —>■ oo limit of the lens elliptic gamma function dill) 


1 — p (pq)^'^+^ 

1 - q-2^ (pq)^^'+i 

“ 1 _ e2i^ q2n (pq)2i+l 

li _ g-2r2 p2n |'pq^2j-|-l 


®It follows from Gi) that a factor TZ{ai, aj, ak) = <l>r,0(1(7? — 2ai))4>7., 0(1(7? — 2aj))4>r, 0(1(7? — 2ak)), would appear 
on the right hand side of the star-triangle relation 1211), if the Boltzmann weights m weren’t normalised by k(q). 


(23) 

n > 0. 
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This function satisfies the following inversion relation 


Q{z,n) = 


1 


Q{-z,-n) 

From this function one defines the edge Boltzmann weights 


(24) 


, X e rnj\ Q(xi — Xj + \a,mi — ruj) Q(xi + Xj + \a,mi + mj) , . 

yVa{cri,aj) = -^^73-7777^- ^ - 73—7 -, (25) 


K{a) 


Q{xi — Xj — ia, rui — mj) Q{xi + Xj — ia, + nij) 


with the normalisation 


K{a) = exp ^ 


-.4072 


n^O 


n((pq)^"' — (pq) 


(26) 


The spectral variable is restricted to the region 0 < a < r], with rj defined in Q. The normalisation 
factor K satisfies the following functional equations 


K{r] — a) 

K(a) 


= Q(i(r/- 2a),0), K(a)K(-a) = 1, 


(27) 


which are required for ([ 6 |) to hold. 

Define also the single-spin Boltzmann weight as 


S{aj) = — Q[2xj — ir/, 2m,) Q{—2xj — iry, —2m,). (28) 

27r 

The continuous spins Xj and discrete spins rrij now take values 

Q < Xj < TT, mj G Z . (29) 

The Boltzmann weights ()25p satisfy spin reflection symmetry 

>V„((Ti, fjj) = W„(cTj, fJi), (30) 


and are 7 r-periodic in the continuous spin Xj. These Boltzmann weights are real and positive for 
p = q*, and 0 < a < rj. 

The Boltzmann weights (125 p satisfy the following boundary conditions 

kVo (<Tj, (Tj ) |^_g — 1, W,^—Q (<Tj, CTj ) |^_^g — r,(p/ r (*5(sin(a;j+ )) rrij T <5(sin(xj J 

(31) 


where S{ai) / 0 . 

The Boltzmann weights (j25p . and (j28p . satisfy the star-triangle relation 

/*7r 

/ dxoS{ao)Wr^-ai{(^i:(^o)^r,-aj{o-j,(To)Wr^-aki(Tk,0-o) = (Uj, (Tfc) (fTj, fJfc) VVa , fJj) , 

moSZ 

(32) 

with r] = ai + aj + ak- 

Note that a similar but different identity, involving an integral and sum over continuous and 
discrete variables respectively, was recently obtained by Gahramanov and Rosengren in the form 
of a pentagon identity from 3 — dJ\f = 2 supersymmetric gauge theories [29]. From this point of 
view, it would also be interesting to find an interpretation of equation (1321) . if any, in terms of a 
new duality for the r —00 reduction of the Si x SsjTjr superconformal indices [^ 1 ^. 
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2.2 Gamma function limit 


From the star-triangle relation (1321) . one can take a further limit to obtain another new solution 
of the star-triangle relation, with Boltzmann weights given in terms of the Euler gamma function. 
The latter star-triangle relation was recently obtained by the author |21j . 

Consider the following limit of the elliptic nomes 

p = e“^, q = e“^, rj = h, h—>0, (33) 


and the following scaling limit of the continuous spins Xj and spectral parameters a from Section 

o 

Xj —>■ hxj, a ^ ha , h —>■ 0 . (34) 

Under the rescaling of the spins aj = {xjh,mj), the asymptotics of (f23]l . in the previous section as 
h —0, are given in terms of the Euler gamma function r( 2 ;) by 


Q{a,) c. {4hr^ 


r( 


l+\raj\+\xj 


r( 


l+\mj\ — \xj 


Then as h —>■ 0, the asymptotics of the Boltzmann weights 
given by 

S{(Tj) ~ ^{Ah)‘^{x] + m^), n{ah) ~ (8h)““ ^ ^ 

ZTT ■' 


(35) 

, and normalisation (121 , are 

1—a \ 


r(i^) 


(36) 


an' 




K{ah)Wahiai,crj) ^ (4fi) 


-4. r(- 


>)r(’ 




l+a+(mi—mj)±\(xi—Xj) ^ l+a+(mi+mj}±\(xi+Xj} ^ 


(37) 

I N 

2 " V 2 / 

In the above limit the star-triangle relation ()32p . formally reduces to the following star-triangle 
relation 

/ CO 

dxo S(ao)Wr,-aii(^i,(^o)yVrj-ajicrj,ao)Wr,-ak((^k,Cro) = Wa ■ (CTj , Ufc) ((Tj, ^fc) Wa^ (fJj, CJi 

mot* 

(38) 

where rj = 1 = ai + aj + anqj 


1 p/'l+o'l -p(i--oi-{mi+mj)±\{xi+Xj)-,-pfl-a-{mi-mj)±\{xi-Xj)-^ 

= _/'™2 , _ 2 ') v\; 12_ 2 __ 2 _; 

' 471-^ ^ ^ P( 1—Q 'j l+a—{mi+mj)±i(xi+Xj) ^ l+a—(mi—mj)±\{xi—Xj) 


-)r(- 


(39) 

The spins of the model now take their values Xj € M, ruj € Z, and the spectral parameter is 
restricted to 0 < a < r], which is a physical regime of the model. These Boltzmann weights also 
obey the spin reflection identity (120^ . however are no longer vr-periodic in the spin. 

The Boltzmann weights p9p also satisfy the following boundary conditions 


^^(crj, iTj)|^_Q — 1, yVri-a{o'iiC!'j)\f^^Q — 2^^^^^{d{Xi-\-Xj) d{Xi Xj) i (4:0) 

“^The following compact notation for products of the gamma function, r(a: ± j/) = r( 2 : -I- j/)r(x — y), is now used 
for convenience. 

typo appeared in the Boltzmann weight S(^aj) in a previous paper | 21| . it has been corrected here by the 
addition of the factor of i. 



















where S{ai) ^ 0. For additional details on this star-triangle relation, the reader is referred to the 
previous publication [2T]. 

Note that different solutions to the star-triangle relation with Boltzmann weights given in terms 
of the Euler gamma function, were previously found in relation to the chiral Potts model |32j . 
The Boltzmann weights for these star-triangle relations are for models that contain only, either 
continuous, or discrete valued spins. It would be interesting to determine if there exists some 
relation between these solutions, and the star-triangle relation given in (I38p . 

As was previously remarked [2T], the appearance of discrete and continuous valued spins here 
resemble the elliptic model obtained by Yamazaki from quiver gauge theory [20] . It has been shown 
here now, how these models are connected through the star-triangle relation (j22|) . the latter relation 
implying the star-star relation for the two-component spin case. 

The alternate method |2T] to obtain (f38ll was to use a scaling limit of the hyperbolic beta 
inteOTal solution of the star-triangle relation [T8| , that resulted in the star-triangle relation given by 
(138p lu The asymptotics of the hyperbolic beta integral in the strong coupling regime, are such that 
sharp delta function shaped peaks appear when the real valued spins take integer values. These 
asymptotics are manifest in the strong coupling limit as additional discrete integer spin variables, 
as appearing in (|38p . One might then ask whether the elliptic variant (|22p of this star-triangle 
relation, arises in the strong coupling limit of some as yet unknown star-triangle relation. Such a 
relation should also then have implications for supersymmetric gauge theories, as well as providing a 
case of an interesting new summation/integration identity, perhaps in terms of more general special 
functions. 

3 Conclusion 

A new solution to the star-triangle relation was given in (|22[) . for an Ising type model whose spins 
contain integer and real valued components. The Boltzmann weights of this model are obtained from 
Yamazaki’s Gauge/YBE correspondence and the related solution to the star-star relation [20|. The 
star-triangle relation (j22p implies the two component spin case of the star-star relation given by Ya¬ 
mazaki. In Appendixj^a proof is presented of a new elliptic hyper geometric summation/integration 
identity which contains the star-triangle relation (j22p as a particular case. The new identity con¬ 
tains six integer variables, in addition to six complex variables, and contains Spiridonov’s elliptic 
beta integral m as a particular case. 

Two further solutions of the star-triangle relation (1321) . (|38p were given that arise as limiting 
cases of (I22p . The Boltzmann weights for these star-triangle relations similarly describe Ising type 
models with integer and real valued spin components. The star-triangle relation (13211 appears to be 
new, while the star-triangle relation (|38h was previously obtained by the author, using a different 
limiting case. In each case the Boltzmann weights are normalised such that Q holds for the 
corresponding lattice model. 

It would be interesting to determine the exact role of the three star-triangle relations (I22p . (j32|l . 
and (l38l) . in the gauge theory setting. The Gauge/YBE duality described by Yamazaki, implies that 
the star-triangle relation (f22]l should correspond to Seiberg duality of the 4 —dA/'=15iX S's/Z^ 
index for SU (2) quiver gauge theory. While the star-triangle relation ([32]) appears to be related to 
the pentagon identity recently obtained by Gahramanov and Rosengren [29], and is thus expected 
to correspond to a duality of indices in 3 — d A/" = 2 supersymmetric gauge theory [30l|3T] . 

It would also be of interest to mathematically prove star-star relations for multi-component 

®Bazhanov, Mangazeev, and Sergeev originally used this limiting procedure to obtain a related solution to the 
star-triangle relation, with Boltzmann weights depending only on the differences of spins. [9] 
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spins given by Yamazaki, that correspond to multivariate generalisations of elliptic hypergeometric 
integral identities of the type in (jA.lSp . It may be possible to do this by adapting proofs given by 
Rains for the continuous variable cases [33], and would likely result in new elliptic hypergeometric 
identities involving an integral and sum over continuous and discrete variables respectively. 
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Appendix A. Proof of ([22]) 


In this section a proof of (1221) is gRen. This is based on and follows closely Spiridonov’s proofs 
of the elliptic beta integrals [141135] kl One major difference is that rather than considering the 
integral over a closed contour encircling the origin, the integral is considered over the interval [ 0 , 27r] 
(the difference between the contours is a simple change of variables). This is done for convenience, 
primarily to avoid calculations involving roots of complex numbers. 

Recall the definition of the elliptic nomes 

p = , q = e™” , = pq , Im((T), Im(r) > 0 , (A.l) 

and now define 

c = i 7 r(l + t/2 - u/ 2 ), (A. 2 ) 

and the following function 

<f{z, m) = {-2ri - 2\z + 2C(|m]l - |-ml)/3) Iml±/(4r), (A.3) 

with z £ C, m r defined in ([ 8 ]), and [mj G { 0 , 1 ,..., r — IT denotes m modulus r|l Define T to 
be the lens elliptic gamma function [ 20 ] in the following forrro 


r(z,m) Y[ 

j,k=0 


I _ e“'^p“l[”^]l(pq)'^+^p”T+i) 1 _ e“'^q“'’+M(pq).?’+iq^(^+i) 
1 - e'^pl[™l(pq)Jp”^ 1 - e'^q”-M(pq)-’q”^ 


(A.4) 


It is useful to introduce the following compact notation for products of this function 


r(x± 2 ;, m±n) := T{x+z, m+n) T{x—z, m—n), r(x±2z, m±2n) := V{x+2z, m+2n) T{x—2z, m—2n). 

(A.5) 

The poles of the lens elliptic gamma function (|A.4I) are located at the points 

z = —TTCJ {rj + |m]) — 2ir]k , —ttt (r(j + 1) — |m]]) — 2\rjk , (A.6) 

^See also Wilf and Zeilberger’s work | 36| . 

®This is equivalent to the definition of \m\r from Section (21 with the r subscript now dropped 

®This is related to 4>(z, m) of dill) by a change of variables 2^; -|- \r], m), and squaring the elliptic nomes 

„ , „2 „ . „2 
P ^ P 1 P q . 
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and its zeros are located at the points 


z = 'Ka {r{j + 1 ) - |ml) + 2\r]{k + 1 ), vrr (rj + |ml) + 2\r]{k + 1 ), 


(A.7) 


for j, A: = 0,1,.... 

The lens elliptic gamma function (IA.4I) obeys the following useful identities 


m) = 


1 


r( 2 i 77 — z, —m) 


and 


n—1 


r( 2 ; + nira, m — n) = 9{z + jvrcr, m — j | r) 1 r(z, m), n = 0 , 1 ,... 

J=o 


where the theta function is defined as 

e{z,m\ r) = q'’)oo(e-'^q 

and 




(A.8) 


(A.9) 


(A.IO) 


4>{z,m) = ^ ^ ^ ^ ^ ^ 2)|m]l± - \{z + 7r)(r - 1 - 2|-m]l)) /(2r). 


(p(z,m) 

The theta function (lA.lOp obeys the following useful identities 


(A.ll) 


0(-z, -m I r) = m | r), 0(z + nrirr, m | r) = ^ | ^ 

(A.12) 

where n G Z. 

A more general identity than the star-triangle relation ()22l) . is the following summation/integration 
identit30 


(q^q'')eo(p^p'')„g d'"^n?=,r(i.±.,.,,±,) 


y=o-^o 


dvr r(±2z,±2y) 


r(tj-1-tj, rtj-b 


u 


•3/ 1 


l<i<j <6 


where 


p,q,tiGC, ttj G Z, |p|,|q|<l, lm(ti)>0, z = l,..., 6 . 

and the variables are restricted to satisfy 

6 6 


^ti = 2i7/, y^ui = o. 


(A.13) 


(A.14) 


(A.15) 


2=1 


2=1 


The star-triangle relation (| 22 p is related to the identity ()A.13p by the change of variables 

h = xi + \ai, ts = X3 + \a3, ts = X 2 - \{ai + - r]), 

t 2 = -xi + \ai, t 4 = -X 3 -|-ia 3 , tg = -X 2 - + a^, - rf), 


and 


ui = mi, U3 = m3, U5 = m2, 
U2 = -mi, U4 = -m3, uq = -m2 


(A.16) 


(A.17) 


° Spiridonov’s elliptic beta integral m corresponds to the r = 1 case of this integral. 
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The identity ()A.13p is what is to be proven. This identity can be re-written in the equivalent 
form 


r2TT 

/(ti,... ... , 7 X 5 ) = V / 

Jo 


dz = 


At: 


y=o -^0 


(q-q’')oo(p^p’')c 


(A.18) 


where 

2/TI) • • • ) ^5) '*^1) • • • ) '^ 5 ) 

and 


, (A.19) 


_ nLi Aiz,Ui± y) r(A -ti,U - Uj) 

T{±2z,±2y)T{A±z,U ±y) ni<j<j <5 '^{U + tj.m + Uj) 

5 5 

A = '^ti, U = '^Ui. (A.20) 


5 5 

= Y^u, 

i=l i=l 

The integral (IA.13h is then recovered by setting tg = 2ir/ — A, and uq = —U. 

The integrand p is 27r-periodic in z 

p{z -P 2Trk, y, ti,..., ts, Ml,..., ns) = p{z, ?/, ti,..., tg, mi, ..., ns), 


(A.21) 


for /c G Z. 

For I Im(A)| < | Im(2ir/)|, the integrand p has the following poles lying in the upper half plane 


{ti + TTU (rj -P Im - 2 / 1 ) -P 2\r]k + 27rn, ti + ttt (r(l -P j) - {m - x/l) -P 2\r]k + 27rn} , 


{-A + vTfj (r(j -I- 1) - |C/ -I- 2/1) -I- 2\p{k -|- 1) -|- 27rn ,-A + vrr {rj + IF" + 2/1) + 2\p{k -|- 1) -|- 27rn} , 

(A.22) 


and the following poles lying in the lower half plane 


{-ti - Tra {rj + fui -P 2 /I) - 2ir/A; -P 27rn, -U - vrr (r(l -P j) - fm + 2/1) - 2ir/A: -P 27rn} , 

[A — ira {r{j -|- 1) — |C/ — 2 /I) — 2\p{k -|- 1) -|- 27rn, A — ttt {rj -|- |f7 — 2 /I) — 2\p{k -|- 1) -|- 27rn} , 

(A.23) 

for i = 1, 2,... , 5, n G Z, and j,k = 0,1, ■■ ■■ By analyticity and periodicity one may also consider 
the integral (lA.lSp (or (IA.13p i. over any contour between the endpoints z = 0,27r in the strip 
0 < Re(z) < 27r, such that the contour separates the points in the two sets of poles (|A.22p and 
()A.23p that lie in this strip, then the ti can be chosen to be any complex numbers as long as such 
a contour exists. 

The idea of the proof is to use a difference equation for p to show that /(ti,..., tg, ni,..., ns) 
is independent of U and n*, z = 1,..., 5, and thus only depends on p and q. Then one can evaluate 
/(ti,..., ts, Ml,..., Ms) using residues at a special value of ti and Ui, to give (jA.lSp . 

The hrst step is to establish the following relation 

I{ti + TT(Tr,t2 ,... ,t5,ni,... jMs) = I{ti ,... ,t5,ni,... ,n5). (A.24) 


To establish this, observe that p satishes the following difference equation 

p{z,y,ti +T:a,t2,...,t:i,ui - I,n2,...,n5) - p{z,y,ti,... ,tii,ui,... ,ut,) 
= G{z - TTa,y + l,ti,... ,ti^,ui,... - G{z,y,ti,. ..,ts,ni,...,ns), 


(A.25) 
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where G is defined as 


G{z,y,ti,... ,t 5 ,ui,... ,U 5 ) = p{z,y,ti ,..., fs,-ui,..., M 5 ) x 

g 2 i 7 rlj/-«i]/r nLi + z,Ui + y\T) Ojii + A,ui + U\t) 

0^2 <^(^1 + Ui+Ui\ t) 0{2z, 2y\T)e{A + z,U + y\T)' 

(A.26) 

To see that (|A.25I) holds, divide both sides of (IA.25P by p{z, y, ti,..., ts, ui,..., us), and one obtains 

e{ti + z, Hi + y I t) e{ti - z,ui - y\T) -p-r 0{A -ti,U -Ui\ r) _ ^ 

0{A + z,U + y\T) e{A - z,U - y\T) e{ti + t*, tti + tti | r) 

^ ei^i/^g(ti + A,m + £/|T) / Q--^^/r+2\Ay-M/r Q{t. + ^ I 

0^2 ^(^1 + k) V 2?/1 t) 6 »(A + z, 17 + y I r) 

iz/r 2i7r(|Tj;-m + ll + |-2?;-ll)/r - Z.Uj - y\T) \ 

e{-2z, -2y \T)e{A- z,U -y\T) j 

(A.27) 

Both sides of this relation are elliptic functions of z, sharing the same poles and corresponding 
residues, then from Liouville’s theorem, the difference of both sides is a constant. The constant 
may be shown to be zero. To check that (IA.27P holds is straightforward, and involves simplifying 

many expressions in terms of |m]]. Some more about the identity ()A.27h is explained in Appendix 

B 

Now integrate both sides of pA.25h . over 0 < 2 < 27r, to obtain 

I{ti + 7rcr,t2, l,U 2 , ...,U5) - I{ti, . . . ,t 5 ,Ui, . . . , 7 X 5 ) 

/ /■27r-i7rIm(CT) (A.28) 

= V / - ]dzG{z,y,ti,...,t 5 ,ui,...,U 5 ), 

y =0 \7-i7rIm((T) JO J 

where 27r-periodicity of G has been used for the first integral on the right hand side. This first 
integral is over a straight line connecting the two points z = —i7rlm(cr) and z = 27r — i7rIm((T). 
Setting Im(tj) > 0, and Im(A) < Im( 7 rr), the function G in (jA.26l) has poles in the upper half plane 
at the points 

{ti + vrcr {rj + [rt* - y]) + 2\pk + 27rn, U + vrr (r(l + j) - fui - yj) + 2\pk + 27rn} , 

{-A + TTcr (r(j + !) + [[[/ + yj) + 2 i? 7 (A; + 1) + 27rn, —A + vrr (r(j + 1) + |t/ + yj) + 2\p{k + 1) + 27m} 

{—A + vrr {r{j + 1) + ft/ + y]) + 27rn, —A + vrr {{U + yl) + 2\T]{k + 1) + 27rn} , 

(A.29) 

and for fU + y} ^ 0 , 

{—A + 7 rr|t/ + yJ + 27rn} , (A.30) 

and the following poles in the lower half plane 

{-ti - vrcj (r(j + 1) + {ui + yJ) - 2iy(/c + 1) + 27rn, -ti - ttt {r{j + 1) - [rii + yl) - 2\p{k + 1) + 27rn} 

{-ti - TTa {r{j + 1 ) + + yl) - 2\r]k + 27m, -U - ttct {rj + {m + y]) - 2 iy(fe + 1 ) + 27rn} , 

{A — 77a {r{j + 1 ) — |t/ — yl) — 2 iy(A: + 1 ) + 27m, A — ttt {rj + |t7 — yl) — 2\p{k + 1) + 27rn} , 

(A.31) 
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and for fui + yj 7 ^ 0 , 


{-ti - nalui + 2/1 + 27rn} , (A.32) 

for i = 1,..., 5, n G Z, and j, A: = 0,1,.... Since for Ini(fi) > 0 and Im(A) < Im(7rT), there are no 
poles in the strip 0 > Im(z) > — Im(7r(T), one can shift the contour of integration in (IA.28P from 
the line Im(z) = — Im(7r(T), to the real axis on [0, 27r], to obtain zero on the right hand side, thus 

I{ti + Tra,t 2 ,... ,t 5 ,ui - 1 ,U 2 ,... ,U 5 ) = I{ti,... ,t 5 ,ui,... ,U 5 ). (A.33) 

Applying this identity r times (as long as parameters remain in our allowed region), one has now 
obtained Equation (|A.24p : 

I{ti + TTar,t 2 , ...,t 5 ,Ui,...,U 5 ) = I{ti,... ,t 5 ,ui,... ,U 5 ). (A.34) 


In a similar fashion, by requiring that Im(A) < Im( 7 r(T), one can establish the following difference 
relation 

I{ti + TTTr,t 2 ,... ,t 5 ,ui,... ,U 5 ) = I{ti,..., Ui, ..., U 5 ) . (A.35) 

For now set Re(cT), Re(T) = 0, Im((T) > Im(r), and rak 7 ^ rrk for any n. A; = 0,1,.... Also let 
the real parts of A, and ti, i = 1,... 5, be non-zero and differ from each other. For the integral 
/(ii,..., As, ui,..., U 5 ) in (lA.lSp . deform the contour of integration from z € [0,27r], such that the 
poles in (]A.22ll . and the following points lie above the contour 


{Ai -|- ITTX 

{—A -|- 2\ri -|- ivrx 

and the poles in ()A.23p . 

{—ti — \ttx I 

{A — 2ir/ — iTTx I 


I X > min (Im(iT)|ui — 2 /l,Im(r)(r — |ni — 2/1)) 

X < max(Im((T)|tti — 2 /l,Im(r)(r — |ni — 2 / 1 )) + 2Im((T)r} , 
I X > min (Im(iT)(r — |[/ -|- 2 /I), Im(r)|[/ -|- 2 /I) — 2\m{a)r A 
X < max(Im(o-)(r -{11 + 2 /I), Im(r)|[/ -b 2 / 1 )} , 
and the following points lie below the contour 
X > min (Im((T) |ui -|- 2/1 > Im(r) (r - |ui -b 2/I )) A 
X < max (Im(cr)|ui -b 2 / 1 , Im(r)(r - \ui + 2/I)) -b 2Im(cr)r} , 
X > min (Im((T)(r — |t/ — 2/I), Im(r)(|C/ — 2 / 1 )) “ 2Im(iT)r A 
X < max (Im((T)(r - {U - 2/I), Im(r)(|/7 - y}))} , 


(A.36) 


(A.37) 


These sets of points correspond to lines in the complex plane with constant real part. Depending 
on the values of Re(Ai) and Re(A), one should translate a set of points by 27rk if needed, so that 
the points always lie in the strip 0 < Re^; < 27r, and the contour of integration remains in this 
strip. Using ()A.24h one then performs n shifts on the variable ti, in the form ti ^ Ai -b vrrr, until 
ti + vrrrn enters the set of points {Ai -b vrarx 11 < x < 2}. Then one transforms Ai ^ ti — nar. 
Under these transformations the poles of the integrand p never cross the contour of integration. 
Thus one obtains 


I{ti + vrrrj - TTark,t 2 ,..., As, ui,..., U 5 ) = /(Ai,..., A 5 , ui,..., U 5 ), (A.38) 

for all j. A; = 0,1,... such that 7 rlm(r)rj — ■K\m{(T)rk G [0, ttI m((T)r]. 
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The set of such points is dense thus I does not depend on ti, and by symmetry on any tj. Then 
from (IA.33I) one has 

I{ti + TTa,t2, . . . ,t 5 ,Ui - 1,U2,. . . ,U 5 ) = I{ti, . . . - I, U 2 , ■ ■ ■ , U 5 ) = I(ti, . . . Ui, . . . , U 5 ) . 

(A.39) 

It follows that I also does not depend on ui and by symmetry on any Uj. Thus I can only depend 
on p and q. 

Set each Ui = 0. In the limit ti +12 ^ 0, p vanishes and the only contribution to the integral 
I, is from two hnite residues coming from poles which cross the contour of integration for y = 0. 
Then evaluating the integral I in this limit from its residues gives the right hand side of (|A.18p . 

By analytic continuation one may then extend the domain of parameter values to that allowed 
by the contour between the endpoints z = 0 , In. 


Appendix B. The identity (IA.27P 

In order to analyse the identity (jA.27p at its poles, it is convenient to write the residues of both 
sides of (|A.27p . in terms of the following theta functioi0 

0(Ar) = (e'^;q^)oo(e-V;q")oo. (B.l) 

Then one hnds that arguments of theta functions appearing on both sides of (jA.27p differ by a 
simple shift Trrrk, A; G Z. The theta function (IB.lh obeys the following useful identity 

6{z + 7rTrk\T) = - — ■ ^ , /c G Z . (B.2) 

^ ^ 0i7rTrfc(fc—1)/2 ^ ^ ^ 

To show that ()A.27p holds requires repeated use of this identity. 

The right hand side of the identity ()A.27p is chosen to have an equal set of poles and residues 
with the left hand side. One wants to show that both sides of (IA.27P define elliptic functions, and 
share the same sets of poles and residues. Then by Liouville’s theorem the difference of both sides 
is a constant, which may be found to be zero. To do this, it should first be shown that no additional 
poles appear on right hand side of (IA.27h . Also it should be shown that both sides are invariant 
under the shift z —)• z + vrrr, and thus define elliptic functions of 2 . 

With the use of (IB.2p . and identities in Appendix C, the calculations involved are straightfor¬ 
ward and are summarised below. 

No additional poles appear on right hand side of (IA.27h ; On the right hand side of (IA.27h . 
there should be no poles appearing at the points 2z = vrr {jr — |—2y]]), or equivalently at 2z = 
ITT {jr + |2y]]), j G Z. By using the identity ()B.2p . at the poles 2z = ttt {jr — |—2y]]) the residues 
of the two terms on the right hand side of ()A.27p will differ by a factor 

exp + ivTufep + rnnkg + AkA + 

for some km,kp,kq,kA, hi G M. It can be shown that this factor is independent of the integer j, and 
is in fact unity, i.e kp,kq,kA-,kti = 0 and km = 0 mod 2. Thus no additional poles appear on the 
right hand side of ()A.27p . 

^^This more resembles the standard theta function appearing in the literature e.g. [,33135] . with a change of variables. 
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Invariance under z ^ z + -KTr: One may use the relation (IB.2P to show that both sides are 
invariant under z ^ z + vrrr. This is the most straightforward property to check. 


Difference of both sides of (jA.27p is zero: Evaluating both sides of (jA.27jl at the point z = 
—ti — tttI—ui — yj, the left hand side gives - 1 , and the right hand side gives 

giil/rgiz/rg2i7r(b-«i + l] + [-2y-l])/r | . _z^Ui-y\T) 

9{-2z, -2y \T)e{A- Z,U -y\T) 0^2 k) 

(B.4) 

After inspection, all theta functions in the above expression cancel up to an overall factor that 
makes the whole expression equal to -1. Thus the difference of both sides of (IA.27P is zero at the 
point z = —ti — 7 rr|—ui — yj. Then since it has previously been shown that both sides are elliptic 
functions of z sharing the same set of poles and conrresponding residues, by Liouville’s theorem the 
difference of both sides is a constant, which must be zero. Thus the identity ()A.27jl holds. 


Appendix C. Useful identities 

Here m G Z, |m] denotes m mod r, and |m]]± denotes mj. 

M± = [-ml± (C.l) 

[—m] + |m — 1] = r — 1 (C.2) 

[m-l]]±-|m]]± + 2[-ml =r-l (C.3) 

[m + 1]± - lmj± + 2[m]] = r - 1 (C.4) 

((2[[m] — r)[m]]± — {2\m — 1] — r)\m — 1]± = —(r — l)(r — 2) — 6[—m] + 6|m]]± (C.5) 

((2[[m] - r)[[m]± - (2[m + 1] - r)|m + = (r - l)(r - 2) + 6[m]l - 6|-ml± (C.6) 
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